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Abstract 

We show higher interior regularity for the Westervelt equation with strong nonlinear damping term of the 
g-Laplace type. Secondly, we investigate an interface coupling problem for these models, which arise, e.g., in 
the context of medical applications of high intensity focused ultrasound in the treatment of kidney stones. 
We show that the solution to the coupled problem exhibits piecewise regularity in space, provided that 
the gradient of the acoustic pressure is essentially bounded in space and time on the whole domain. This 
result is of importance in numerical approximations of the present problem, as well as in gradient based 
algorithms for finding the optimal shape of the focusing acoustic lens in lithotripsy. 
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1. Introduction 

High intensity focused ultrasound has numerous applications starting from the treatment of kidney and 
bladder stones, via thermo therapy and ultrasound cleaning to sonochemistry. Due to the nonlinear effects 
observed in the propagation of ultrasound in these cases, such as the appearance of sawtooth solutions, 
models of nonlinear acoustics and their rigourous mathematical treatment have become of great interest in 
recent years. 

One of the most popular models for the nonlinear propagation of ultrasound is the Westervelt equation 

(1 — 2fcu)ii — c^Au — 6Au = 2fc(u)^, (1.1) 

expressed here in terms of the acoustic pressure u, where a dot denotes time differentiation, b the diffusivity 
and c the speed of sound, k = /3a/A, A = gc^ is the bulk modulus, g is the mass density, /3a = 1 + B/(2A) 
and BJA represents the parameter of nonlinearity. A detailed derivation of (11.11) can be found in [l0|, [ISj 
and [ 2 ^. 

Westervelt’s equation is a quasilinear wave equation which can degenerate due to the factor 1 — 2ku. 
This means that any analysis of this equation has to include bounding away from zero this term, i.e. finding 
an essential bound for u. That has been so far achieved by means of employing the Sobolev embedding 
H^{Q) ^ (cf. [Ilj, [13)) which implies that the solution of the Westervelt equation has to exhibit 
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regularity in space. However, achieving TJ^-regularity is too high of a demand in the case of coupling 
acoustic regions with different material parameters. 

To remedy this issue, Westervelt’s equation is considered with an added nonlinear damping term 


(1 — — c^Au — div(6((l — 5) + ^1 ^)Vu) = 2fc(ft)^, (1-2) 

cf., i, [13, where <5 S (0, 1), q > 1, g > d— 1, and d G {1, 2,3} is the dimension of the spatial domain H C 
on which (HU is considered. Since degeneracy is avoided with the help of the embedding ^ 

L°°{n), using this model allows to show existence of weak solutions with regularity in space, and in 

turn well-posedness of the acoustic-acoustic coupling problem. 

Considerations of the acoustic-acoustic coupling are motivated by lithotripsy where a silicone acoustic 
lens focuses the ultrasound traveling through a nonlinearly acoustic fluid to a kidney stone (see fig. 1). 
The interface coupling is modeled by the presence of spatially varying coefficients in the weak form of the 
equation (II.2p (see [3 for the linear and [3 and 171 for the nonlinear case) as follows: 


Find u such that 

Iq - ‘^Kx)u)u4> -b ■ \/(j) + b{x){l - 6{x))\/u ■ Vcj) 

+b{x)S{x)\Vu\‘^~^Vu ■ V4> — dxds = 0 

holds for all test functions (j) € X = L^(0,T; lFg^’‘^^^(r2)), 


with (tt,'d)|t=o = (wo,ui). In this model b stands for the quotient between the diffusivity and the bulk 
modulus, while the other coefficients maintain their meaning. The coefficients are allowed to jump only over 
the interface, i.e. the boundary of the lens. For notational brevity, we emphasized the space dependence of 
coefficients in (ra, while omitting space and time dependence of u in the notation. 


Axis of Rotation 

Kidney stone 



fig. 1: Schematic of a power source in lithotripsy 
based on the electromagnetic principle 


The first goal of the present paper is to show higher interior regularity results for solutions of lO]). 
We will show that u G and |Vm|^Vm G Althou^ q-Laplace and 

parabolic g-Laplace equation have been extensively studied in the past (see [13 , [13 , 0 , 13 references 
given therein), regularity results in literature on hyperbolic equations with damping of the q-La plac e type 
are sparse and have so far been concerned with local and global well-posedness (see [13,0, i,113)- 
Secondly, we will consider the coupled problem and show that the solution to m is piecewise regular 
in space under the assumption that the gradient of the acoustic pressure remains essentially bounded in space 
and time. This result is crucial in future numerical approximations of the present problem, as well as in 
gradient based algorithms for finding the optimal shape of the focusing acoustic lens, where regularity 
of u is needed in order to express the shape derivative in terms of integrals over the boundary of the lens, 
see [ 13 . 
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1.1. Overview 

The paper is organized as follows. We begin in Section[2]by recalling certain basic results from the theory 
of finite differences, as well as some helpful inequalities which are needed later. In Section[3]we prove higher 
interior regularity for the Westervelt equation (11.21) . Section 0] extends these results to the coupled problem. 
In Section [SI we will show that we can obtain higher regularity up to the boundary of the subdomains if the 
gradient of the acoustic pressure is essentially bounded in space and time and the subdomains are sufficiently 
regular. 


2. Preliminaries 

In what is to follow, we will need to employ difference quotient approximations to weak derivatives. 
Assume that C K.'^, d G {1, 2, 3} is an open, connected set with Lipschitz boundary. Let V CC fl. will 
stand for the r-th difference quotient of size I 

for a; G fo, Z G R, 0 < |Z| < idist(fo, dfl). Then D^u := {D\u ,..., D^u). We recall the integration by parts 
formula for difference quotients 

/ uD^j.ipdx = — / D~’‘ufdx, 

Jv Jv 

where if G C^{V), 0 < |Z| < ^dist(fo, 9S1), as well as the product rule 

Dl{fu) = f^Dl.u + uD^f, 

with f^{x,t) ■= fix + Icryt). We will also need the following result (cf. Theorem 3, Section 5, Q): 
Lemma 1. (a) Assume 0 < q < oo and u G IT^’'^'''^(r2) Then for each V CC Tt 

||L>*u||/,9+l(y) < CllVwllig+qo), 


for some constant C and all 0 < \l\ < ^distiV, ^^l). 

(b) Assume 0 < q < oo, u G L^'^^iCl), and there exists a constant C such that ||T*^n||L<z+i(y) < C for all 
0 < |Z| < ^dist{y,dTl). Then 

u G with llVwllig+qy) < C. 


Essential inequalities 

Before proceeding further, let us also recall several useful inequalities that we will need when handling 


the (j-Laplace damping term. They can be found in Chapter 10, [l^ and Appendix, [l^. From now on, C, 


will be used to denote a generic constant depending only on q. For any x, ?/ G R“ it holds 


I"* x-lvl"^ y\ < Cg\x - y\i\x\ + \y\y , q > 0, 


( 2 . 1 ) 



\\xr^x-\y\<>-^y\>2-^\x-yW\x\ + \y\Y-^ 

> 2^"«|a: - y|«+^ >0, y > 1. 

(2.2) 


(q + 1)2 11^1 " 2 ; \y\ ^ yp < (|a;|« ^x 

1 

1 

y), <? > 1, 

(2.3) 


\\x\'^-^x - < <?(|a;|'^ + \y\'~^) 


, <? > 1- 

(2.4) 

We will also 

need Young’s inequality (see for instance Appendix B, Q) in the form 



|a;y| < e|a;r + (^(e, r)|y|--i 

(£ > 0, 1 < r < 

oo). 

(2.5) 

with C{e, r) 

V 1 

= (r - l)r^£~^. 
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Notation 

By Cx Y denote the norm of the embedding operator X (O) —?> Y (O) between two function spaces 

over the domain fl. 

3. Interior regularity for the Westervelt equation with strong nonlinear damping 

In this section, we will establish higher interior regularity for the equation (1121) with constant coeffcients. 
Let us hrst consider the following Dirichlet problem: 

{ (1 — 2ku)u — c^Au — 6div(((1 — d) + <5 |Vm|‘^“^)Vm) = 2k{u)^ in x (0,T], 
n|an = 0 fortG(0,r], (3.1) 

(M,u)|t=o = (uo,ui) in O, 

with the following assumptions on the coefficients and the exponent q: 

c^, b > 0 , (5 G (0,1), k gR, q > d — 1, q > 1. (3.2) 

The weak formulation reads as 

{ /o^ /o{(l - 2ku)u(j) + c^Vm • Xcj) + 6(1 - (5)Vm • Xcj) 

+b6\Vu\‘>-^Xu ■ V(t> - 2k{uf(t>} dxds = 0 (3.3) 

holds for all test functions (j) G X = L^(0, T; Wg^’^''"^(n)), 

with (u,m) = (uo,wi). We recall the following well-posedness result (cf. Theorem 2.3, i): 

Proposition 1. (Local well-posedness) Let assumptions (13.21) hold. For any T > 0 there is a kt > 0 such 
that for all uo,ui G with 

l^ilL2(n) + l^“olL2(f2) + |Vui|^2(Q) + + l^“oli9+i(n) ^ 

there exists a weak solution u G W C X of (13.111 . where X = iL^(0, T; L^(r2)) (~l (^^’^(O, T; Wg^’‘^~'’^(n)), and 

W = {v G X : ||h||i2(o,T;L2(f2)) ^ ^II^^IU“( 0.T;L2(n)) < ™ 

A||Vb||L,+i(o,T;Li+i(n)) < M A{v,v) = (uo,ui)} (3.4) 

with 

2\k\C^,,,+,^^^{KT + T^M) <1 (3.5) 

and m sufficiently small, and u is unique in W. 

In 0, the issue of possible degeneracy of the Westervelt equation due to the factor 1 — 2fcM is resolved 
by means of the embedding ^ L°°(Ll), valid for q > d — 1, and the following estimate 

\u{x, t)| < C'^c9+i_ioo|VM(t)|Li+i(n) 

-^wi.9+1 looIVmq + f VMds|i9+i(o) 

“ ’ Jo 

pi p 

<C'|^c<.+i_ioo(|VMo|L9+i(n) + jjS/u(y,s)(‘+'^dyds^‘‘*^^, 

which leads to the bound 

1 — ao < 1 — 2ku < 1 + oo, 

/O f{\ 

ao := 21^1(7^1,(iV^olLi+qn) + Ti+^ || V'u||i;,g+i(o,T;L9+i(n)))- 

Due to the embedding lT^’‘^+^(n) ^ C^J~1TT also know that u is Holder continuous in space, i.e. 

u G COA(0,T; C'°’^”i^(H)). 
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3.1. Higher interior regularity 

We will establish higher interior regularity by following the difference-quotient approach (see, for instance, 
Theorem 4, Section 6.3.2, Q). Let us denote 

F=|Vu|^Vm, = IVm'I^Vu'. 

As a by-product of the following proof we will also obtain F G by adapting the idea of Bojarski and 

Iwaniec for g-harmonic functions (see, for instance, Section 4, |15| ') to our model. 

Theorem 1. (Higher interior regularity) Let assumptions (13.21) hold true, ug G H^(fl) D Wg’'^'^^{V l), 
Ml G Wg^’^~'"^(n), and let u be the weak solution of (13.111 . Then u G H^[0,T■ Hf^^{n)) and |Vm|V^Vm G 
L^iO,T;HUn)). 

Proof. Choose any open set V CC fl and an open set W such that V CC W CC Ll. We then introduce a 
smooth cut-off function f such that 


I C = 1 on F, ^ = 0 on 11 \ W, 

jo<c<i. 

Let 1^1 > 0 be small and choose r G {1,..., d}. We are then allowed to use 

0:=-Zl-'(C'4M)X[o,t), tG[0,r] 

as a test function in (13.31) . which results in 




(1 — 2kw){fD(uY dx 


o”^ 2^ 


|CL>J.VmP dx 




J 0 


+ b{l-6) [ [ \(:Dl.S/u\^dxds + -^^b6 [ [ \CDlF\^dxds 

Jo Jdi \Q + 1)'^ do Jn 

<2k f f D^uuCYD^iidxds — (? f f CfVC, ■ u D\u dx ds 

Jo Jo, Jo Jo, 


-2b{l-S) 


CVf ■ D)\/u D^iii dx ds 


(3.7) 


0 Jo 


— 2bdi f f Z1(,(|Vm|'^ ■ (VC Dliidx ds + k f f [if + 2u){CD^^uY dx ds. 

Jo Jo Jo Jo, 


/o Jo 

Here we have made use of the estimate 
ct 


CD){\Vu\<i-'^Vu)D)Vudxds > 


0 Jo 


(9 + 1)^ 


1 . 2 11Y7 ■ l \ — ■ ’ '— . ' . 


0 Jo 


p 


Cl|Vu'| —Vm‘ - \Vu\ — Vu\^dxds, 


which follows from dZSl). Next, we estimate the terms on the right hand side containing C^C- We have 


— c 


[ f CVC-DlVuDl.udxds-2b{l-6) [ [ (Vf ■ DlVuD^iidxds 
Jo Jo Jo Jo 


<C 


C{\DlVu\ + \DlVu\)\VCDlu\dxds 


10 Jo 


C. 


< er||CLlrVM||^oo(o_T;L2(0)) + £||CA>rVM||j;^2(o^T;L2(n)) + — II^^IIl2(0,T;L2(O)), 

where C depends on c,b,5 and |VC|l=o(w) and we have used Lemma [I] (a). By employing estimate (j2.4l) 
and Holder’s inequality we obtain 


-2b6 [ f Dl{\Vu\'^-^Vu)-CVCDludxds 
Jo Jo 
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< 


2b5 [ [ y|D[u| - |Vu|«-^Vw|C|VC|(ia:ds 

Jo Jn ‘ 

rT 


< 2qb5 


10 Jn 

rT 


|VC£>[u|(|Vm'|^ + iViil'^^lCDlFldxds 


<2qb6( f f \VCDI.u\'‘~^^ dx ds\'’^^ ^ f f dVu^l’^ +|Vu|”*2 ) ' dx dsl 

^Jo Jn ^ ^Jo ^suppc ^ 


q-i 

2(<I+1) 


Uo I 


10 Jn 

The second integral can be majorized with the help of Minkowski’s inequality by 


{/ f +|Vm|'*2 ) 9-1 'dxdsj 

^Jo JsuDvC ^ 

dxds^ 

^lo =2||V^ 


g -1 

2(g+l) 


0 «/suppC 

T 

0 «^suppC 
T 


g-i 

2(g+l) 


1 / f {\^dx ds\ 

^Jo JsuddC ^ 


g-i 

2(g+l) 


0 «/suppC 

g-i 

2 

_Lg+i(0,T;L‘2+i(n))’ 


for small |/|. Utilizing Young’s inequality then yields 

-2bS [ f Dl{\Vu\'i-^Vu) ■ (VCDliidxds 
Jo Jni 

<b5e j^\(:DiF\^ dxds + 

where C > 0 depends on bi, 6i, e, q and |VC|l~(vk)- Note that the first term in the last line can be absorbed 
by the biSi- term on the left hand side in (ITtI) for sufficiently small £ > 0. The two remaining terms on the 
right hand side in (EJl) can be estimated as follows 

k [ [ (u^ + 2 u)(^{DIu)^ dx ds + 2k f f u D\u D\u dx ds 
Jo Jn Jo Jn 

< C'(ll'“llL“(0,T;L°°(a))IIC-DrM|li2(o,T;L2(M^)) 

+ (C'Si,L4)^lkllL2(O.T;L2(n))||C-DrU||ioo(o,T;ffi(n))||C^r^||L2(o,T;£ri(n))) 
^C'(ll^lU“(0,T;L°°(n)) \\Dr'k\\L^{o,T\L^{n)) 

+ ™^^(ll'Or'“lli“( 0 ,T;L 2 (vi/)) + II II i“( 0 ,T;L 2 (O)) ) 

+ £||T)J,m|||2(o_T;L2(vI/)) + £|IC-DrVM|||2(o_T;L2(n))) • 

Altogether, for sufficiently small £ > 0 and m, we can achieve that 

IIC-C*r^llL=°( 0 ,T;L 2 (O)) + IlC-C^r^^i IIl=°( 0 ,T;L 2 (O.)) + II ^ IIL 2 (q_- 2 -. 2,2 (q)) 

+ IIC-Dr-F|li2(o_T;L2(o)) 

<C(||Vhr^t\o,T;L~(n)) + P'«ll L9 +i(0,T;L9+i(W)) II^IU“(0,T;L“(n))||-Dj,u|||2(o_T;L2(w)) 

+ ll'^r'“lll,“(0.T;L2(W)) + ll'^r^llL2(o,T;L2 (w)) + l'^r^llL2(0) + |^2 (f;) ) . 

By remembering the definition of C and Lemma [U we finally arrive at 

||-Drni|li°°(0,T;L2(y)) + l|-C>[Vu|||oo(o^T;L2(y)) + l|-C>[V'u|||2(o^T;L2(y)) + ll-^r^lli2(0,T;L2(y)) 









< C(||VM||^'to\o,T;L“(n)) + (1 + ll«llL“(0,T;L“(n)))ll Vu||i2(o_T;L2(n)) + ll^^^llL“(0,T;L2(n)) 

+ I^^ili2(n) + l''^o|ff2(Q)), 

for r e [1, d], sufficiently small |?| > 0 and sufficiently large C > 0 which does not depend on Z. By employing 
Lemmadl we can conclude that u £ H^{0,T-, and |Vu|^Vm G L^(0, T; □ 

As a simple consequence of the previous proposition, we can obtain Holder continuity of u (see Section 
4, [l^). Indeed, since F G L^(0, T; and d G {1, 2, 3}, due to Sobolev’s embedding theorem we have 

that F G L^(0, T; Lf^^(r2)). This implies that u G We can than conclude that 

u G Wi-^+i(0,T;C'°’“(H)), where a = 1 - 

When d G {1, 2} we can do even better. According to Sobolev’s embedding theorem, u G (0, T ; (H)) 

if d = 1, and u G if d = 2, where 7 G (0,1). Altogether, we have 

Corollary 1. Let the assumptions of Theorem^hold true. Then 

ifd = 3, 

uG I Wi-9+i(0,r;C°’J(H)) mtd7€(0,l), */d = 2, (3.8) 

[t4^M+i(0,r;C;i(H)) ifd=l. 

3.2. Neumann problem for the Westervelt equation 

Let us also consider the Neumann problem for the Westervelt equation with strong nonlinear damping: 

(1 — 2ku)u — c^Au — 6 div(((l — d) + d|V'u|‘?“^)Vu) = 2k{u)^ in H x ( 0 ,T], 

c2f^+5((l-d) + d|Vu|«-i)f^=5 on5Hx(0,T], (3.9) 

{u,u)\t=o = iuo,ui) onH,, 

with the same assumptions (13.21) on coefficients. Problem (13.911 is locally well-posed thanks to the following 
result (cf. Theorem 2.5, 0 ): 

Proposition 2. Letg G L°°{Q,T\ {9^))) 9 S L^ (0,T; (i^f^)): anduQ.ui G 

For sufficiently small initial and boundary data, final time T and m and M there exists a unique weak so- 

lutionuGWcX of dSH), w/iere a: = ij2(0,T;L2(H)) nC°d(0,T;ITi'«+i(H)), and 

W = {v £ X : ||h|lL 2 (o,T;L 2 (o)) < m A ||u||icx=(o,T;Hi(n)) < rh 

A IIVh||i9+i(o,T;L9+i(f2)) — ^}- 

By inspecting the proof of Proposition [U we immediately obtain higher interior regularity result for the 
present model, since the cut-off function used in the proof vanishes near the boundary: 

Corollary 2. Let the assumptions (13.21) hold true, uq G H^{Q) l~l ui £ and let u 

be the weak solution of (13.9p . Then u £ Lf^(0,T; and iVul^Vu G L^(0, T; iLj^^^(H)). Moreover, 

(I3.8|l holds. 

4. Interior regularity for the coupled problem 

Let us now assume that ft C M'’*, d G {1,2,3}, is a bounded domain with Lipschitz boundary dfl, and 
a subdomain, representing the lens, such that 17+ C Ll and H+ has Lipschitz boundary cA7+ = P. 
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We denote by H.- = il\ f2_|. the part of the domain representing 
the fluid region. We then have = F U dQ. 
n_|_, n_ will stand for the unit outer normals to lens and 
fluid region f2_. Restrictions of a function v to will be 

denoted by u+, V- and |z;] := v+ — V- will denote the jump 
over r. 

Note that the assumption on the regularity of the subdomains 
will be strengthened to when showing higher regularity up 
to the boundary of the subdomains. 

Lens and fluid Q— regions 



(4.1) 


The coeflicients in dLSD will only be allowed to jump over the interface F: 
f b, Q, X,S,k e 

\bi := b\Q.,gi := := > 0, := £ (0,1), ki := k\n. £ R for i £ {+, -}. 

We assume that q > 1, q > d — 1. The strong formulation of reads as follows: 

^(1 - 2k{x)u)u - div(^Vu) - div^&(a:)((l - (5(a:)) + (5(a:)|Vu|‘?“^)Vu 
= inn+Ufl_, 

|m] = 0 on F = 9n+, 

I Jit + =0 on r = 

u = 0 on 

.(■u,u)|t=o = (wo,Mi)- 

This model was studied (in an equivalent one domain formulation) in [^. We will utilize the following 
well-posedness result (cf. Theorem 2.3 and Corollary 4.1, i): 

Proposition 3. (Local well-posedness) Let q> d—1, q>l and the assumptions (14.1|) hold. For any T > 0 
there is a kt > 0 such that for all uq, ui £ WQ^’'^~'"^(f2) with 


(4.2) 


I^ili2(n) + I^^^oli2(o) + |VMi||2(f2) + |Vuoliq+i(n) + < n 


9+1 < ^2 

(n) - 


there exists a unique solution u € W C X of (14.21) . where X = iF^(0, T; L^(r2)) D(^^^(O, T; WQ^’'*~'’^(n)), and 


W = {vGX 


with 


l|ii||L 2 (o.T;L 2 (o)) < mA ||V'y||i=o(o_T;L 2 (n)) < m 
A||Vf;||i,+i(o,T;Li+i(n)) < M A{v,v) = (uo,wi)}, 


2kC^,,,+,^^^{KT + T^^M) < 1 , 


(4.3) 


and m sufficiently small. 


For simplicity of exposition, higher interior and later boundary regularity will be obtained under the as¬ 
sumption that the coefficients in (14.21) are piecewise costant functions, i.e. 


h := b\n ^, pi := , Ai := := , h := k\n^ are constants, 

bi,gi,Xi>0, (5i£(0, 1 ), fci £ K for i £{-(-,-}. 


(4.4) 


We denote a; = min{|cLi+|, |a;_ |}, uJ = max{|a;+|, jw-1}, where w £ {&, p. A, <5, k}. The proof of Theorem [T] 
can be carried over in a straightforward manner to the coupled problem to show higher interior regularity 
within each of the subdomains: 











Corollary 3. Assume that q > I, q > d — I, uo\n^ € i G {+, —}, uo,ui € (fl) and assump¬ 

tions (03) on the coefficients hold true. Let u be the weak solution of (Ol) . Then Ui e H\Q,T-Hl^{Vti)) 
and G L^{0,T; i G {+,—}. Moreover, (13.81) holds with u replaced by Ui and Ll by 

Lli, i G {+, —}. 

4.T Neumann problem for the coupled system 

Let us also consider the coupled problem with Neumann boundary conditions on the outer boundary of 
the fluid subdomain: 


- 2k{x)u)u - div(^VM) - div(6(a:)((l - S{x)) + (5(x)|V'u|‘? ^)\7u) 

= inf7+Uf7_, 




[u] = 0 on r = dfl+, 

'^iffi + ba-S)^ + bS\Vu\ 


a—1 du 

drij^ 

— 1 du 


= 0 


dn 


= g 


on r = 
on dTl, 


(M,u)|t=0 = (mo,Mi)- 


(4.5) 


We will show that the higher regularity result is valid for this model as well. The weak form of the problem 
is given as follows: 



i(l - 2ku)u(t) + -Vu ■ V(/) + 6((1 - ,5) + (5|Vm|«“1)Vu • Vtj) 
X g 




dx ds 


for all (j) G L^(0, T; with initial conditions {uq, ui). We recall the following well-posedness result 

(cf. my- 

Proposition 4. Let assumptions f|4.ip hold. Let q>d — 1, q>l, gG L°^{0,T;W~^'~^ {dLt)), g G 

L~ (dfl)), and uq,ui G For sufficiently small initial and boundary data and 

final time T, there exists a unique weak solution u G Vd G X of (|4.5I) . where X = iJ^(0,T; L^(r2)) D 
and 


W — {v G X : ||'c||l 2 (o_T;L 2 (o)) < m a ||f'||L“(0,T;Hi(f2)) - ™ g\ 

^ IIV^llL9+i(0,T;L9+i(n)) < M}. 

The interior regularity result can again be transferred from Corollary [3] to the present model: 

Corollary 4. Let the assumptions of Proposition 0] and assumptions (14.411 on the coefficients hold true, 
let Moloi G i € {+!“}; uo,ui G and let u be the weak solution of (14.5|) . Then 

Ui G H^{0,T-, and \Vui\^^Xui G L‘^{Q,T-,Hl^Jffili)), i G {+,—}. Furthermore, (13.811 holds with u 

replaced by Ui and ft by fti, i G {+, —}. 


5. Boundary regularity for the coupled problem 

We will show next that the iJ^-regularity result can be extended up to the boundary of each of the 
subdomains under the assumption that the gradient of u is essentially bounded in time and space on the 
whole domain. For this property to hold, we will need to smoothen out the subdomains, i.e. assume that 
they are regular. The proof will expand on the approach taken in Lemma 3.6, Q. 
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Theorem 2. (Boundary H'^-regularity) Let the assumptions of Corollary^ hold and let dfl and F = dfl 
be regular. If u G T; VF^’“(f2)) and ||VM||i °°(n)(o T-L°°(n)) sufficiently small, then u+ _ 

Proof. We will only show that u+ G H^{0,T; H^{Q+)), since U- G iF^(0, T; iJ^(n_)) follows analogously. 

Step 1: Straightening the boundary. We begin by straighteniM the boundary through the change of coordi¬ 
nates near a boundary point (cf. Theorem 4, Section 6.3.2, [8|). Choose any point xq G dLl+. There exists 
a ball B = Brixo) for some r > 0 and a C^’^-diffeomorphism ; i? —>• C such that det|V4'| = 1, 

U' = '^{B) is an open set, n f2+) C Mlf. and n F) C where R-lj. is the half-space in the new 

coordinates. 


$ 



We change the variables and write 

y = 4'(a;), x G B, 

X = $(?/), y G U'. 

Then we have 4'(i? n f2+) = {y € C/' : j/„ > 0}. We denote 

B+^Brn n+, G = 4/(5. (xo)), G+ = 4/(5+). 

Then G CC 5' and G+ C G. We define 

w{y,t) := u{^{y),t), {y,t) G U' x [0,T]. 

It immediately follows that w{t) := w{-,t) G We now transform the original equation on 

5 X [0,T] into an equation on U' x [0,T]: 


1 ^ 

l-{l-2kw{t))w(t)(j)+ ^ (^^cFijDiw{t)Dj(j) + ^ijDiWi{t)Dj4> 


lu' '' X 


ij=l 


lJlVw(t)l'^ ^f]ijDiw{t)Dj(t^ - -j-{w{t))'^(l)j dy = 0, 


for a.e. t G [0,T], and all (f G where DiW = and 


(5.1) 


X{y) = A($(y)), k{y) = k{<I>{y)), 

d 

a. 






^ p(4>(y)) dxr dxr 

d 


4(y) = 

r=l 


(5.2) 
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(Ti + 





r=l ^ ^ 

Note that Draij,Dr^ij,Drfjij € L°°{U') for r S {!,...,d — 1} since and $ are mappings and 
= ] 4 / 2 ,oo Chapter 2, Section 2.6.4, Q). It can be shown (cf. Section 6.3.2, Q) that 

d d 

'y ^ ^l|^l ) ^ ^ 5 

i,i—l ^,7 = 1 

(5.3) 

> Ki\cj)\'^, gU' xR'^. 

i,j=i 

Next, we choose a domain W' such that G CC W' CC U' and select a cut-off function such that 

f C = 1 on G, C = 0 on \ II^', 

jo<c<i. 

Let |Z| > 0 be small and choose r G {1,..., d — 1}. Note that since we consider directions parallel to the 
interface now, we have D^X = 0, = 0 and that there exists a constant K 2 > 0 such that 


\Dla,j{y)\ < K 2 , \D%{y)\ < ^ 2 , \D%j{y)\ < (5.4) 

for a.e. y G W, I <i,j < d and sufficiently small |Z|. 

Step 2: Existence of second order derivatives DjDiW G L'^{G'^)), j 7 ^ d. We then use (j) = 

—D~^{Cf^D\.w{f)) as a test function in (15.11) . which, after integration with respect to time, results in 


1 r /■ 1 


1 ‘ 1 


2 ^Jw A* 

d 


{l-2k''w'^){C,Dlwf dy D[DiwD'-^DjW dy 


0 2 Ljvf' ““1 

%,3 — l 


-1 
- 0 


f f ^ nlDiwDlDjW dy ds 

Jo 


‘ 2 /c' , 


ft f fl 

dyds+ / —{w’'+2w){fDlw)^dyds 
Jo Jw A* 


Dfw ( 

0 JW’ A^ 

pt p ^ 

/ / Dl{aij)DiW {(f^D^DjW + 2 CDj(DIw) dyds 

Jo 


2 j f aljDlDiW (DjfDlw dy ds 

Jo 

pt p ^ 

/ / y^ Dl{iij)DiW (f^DlDjW + 2fDj(Dl.w) dyds 

Jo 

ft n d 

2 / / y^ fijDiW CDjfDlw dy ds 

do J W' 1 


/o 

d 


[ f ^ DI.{\JIVw\'^ ^fiijDiw){C‘^Dl.DjW+ 2CDjCDl.w)dyds. 
Jo 
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We can estimate the last term on the right hand side as follows 


pt p ^ 

- Dl.{\JlVw\‘^~^fiijDiw){('^Dl.DjW + 2({DjC)DIw) dyds 

Jo Jw'i j^^ 

Jo 

pt p ^ 

- Dl{\JlVw\‘^~'^fiij)Diw{C'^DlDjW + 2({Dj()Dlw) dyds 
Jo Jw ■ ■ 1 

d 

<C||Vzi;||r4,T;L~(w'))(EllC^'^*^lli=(0, T-,L^{W')) + \\'^C\\l°°{0,T;L°°{W'))\\^w\\'l'^(0,TiL^{W'))) 


with C independent of VC- Due to (15.41) the rest of the terms on the right hand side can be estimated 
analogously to the estimates in the proof of Theorem [U which for sufficiently small ||Vr(;||ioo(Q_ 7 -.j;^oo(i 4 //)) 
leads to 


d d 

{0,T-,L^{W')) +E IICD[Di'U;|||oo(0^T;L2(n^/)) + \\C^r^iM\'L2{0,T-L'^{W')) 

i=l i=l 

<C'((1 + II^IU“(0,r;L“(lV')))ll^^lli“(0,r;L2(iy/)) 

+ (1 + II ooj-Q Vffi||j;^ 2 ('Q y 

■L'^(W')) + |Vu;(0)|^2(vi//) + |ii^(0)|^2(vF'))' 

Recalling the definition of C and employing Lemma [T] yields DjDiW € iL^(0, T; L^(G'+)) for 1 < z < d, 

1 < j < rf - 1 - 


Step 3: Existence of second order derivative DdDdW G JI^(0,T; L^(G~^)). It remains to show that DddW := 
DdDdW G id^(0,T;L^(G+)). From (15.11) . after integration by parts, we obtain 



{uddDdw{t) + iddDdwit) + Ij|’Vu;(t)|‘? ^f]ddDdw(t)]Dd(i)dy 


IJ 


1 - 9 k 

- {l- 2 kw{t))w{t) + ^{w{t)r 

A A 


d— 1 d 

+ '^'^(^Dj{a^jD,w{t)) + Dj{iijDiw{t)) + D^dj|’Vw(t)|'^“^%Affi(O)) 

1=1 i=l 


(5.5) 

(5.6) 



f{w){t)(j)dy, 


for (j) G G“(G+), a.e. in [0,T]. Since the right hand side of the equation is well-defined, we conclude that 
for a.e. t G [0,T] the weak derivative of addDdw{t) + fddDdw{t) + \J^S/w\'^~^fiddDdw{t) with respect to yd 
exists on G"*". Furthermore, for a.e. t G [0,r] the weak derivative satisfies 

- Dd(crddDdw{t) + iddDdw(t) + \JlS/w\'^~'^f}ddDdw{t)) = f{w){t) (5.7) 


on G"*". From what we have shown, it follows that f{w) G L^(0, T; L^(G''')). We set 

z{t) := addDdw{t) + iddDdw{t) + \JlVwf^~'^'qddDdw{t), 


and 

fdd{t, Ddw{t)) := fdd + |d$ Vw(t)|'^“^77dd, 
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(suppressing in the notation dependence on Diw{t),... ,Dd-iw{t) which we already know to be smooth 
anyway) so that relation (15.71) reads 

-Ddz{t) = f{w){t), 

where 

zit) = addDdw{t) + ^ddit, Ddw{t))Ddw{t). (5.8) 

Since f{w) G L^(0, T; L^(G+)), and using the fact that DjDiW G iJ^(0, T; L^(G+)) for 1 < i < d, 1 < j < 
d— 1, we have 0 G On the other hand, due to assuming that u G W^^’°°(0, T; fF^’°°(r2)), 

we know also that z G L°“(0, T; L^{G~^)). Therefore we conclude that 

z G L^{0,T;H\G+)) nL°°iO,T; L“(G+)). 


Since dSH) represents an ODE (pointwise a.e. in space) for Ddw{t), it can be resolved as follows: 


Ddw{t) = exp^- 


<^dd 


■ dr 


z(r) 


■ exp 


CTdd 


0 ^dd{T, Ddw{T)) / Vdo ^dd{T,Ddw{T)) Vdo ^dd{p,Ddw{p)) 
+ Ddw{Q)^, 

and then also Ddw{t) can be expressed in terms of z,add,^dd' 

zit) ddd 


dp'j 


dr 


(5.9) 


Ddwit) = 


^ddit,Ddwit)) ^dd{t,Ddw{t)) 
z(r) 


■ exp 


0 UdiT,DdwiT)) \Jo ^ddip,Ddwip)) 


exp - 


<^dd 


f^dd 


0 ^ddir, Ddwir)) 


■ dr 


dp I dr + DdwiO) 


(5.10) 


Since the right hand side depends on Ddwit), this is rather a fixed point equation than an explicit expression 
for Ddw{t). We therefore consider the fixed point operator T : M —)• M, defined by the right hand side of 
(I5T01) . i.e., 

T = Son, 

where TZ : M ^ M, M = {(p G M : (p > Ki}, is the superposition operator associated with ^dd, i-e., 
TZiv)it) = ^ddit,vit)), and S : M ^ M 


c/ ^i^) ddd f f ddd 

‘5(r)(t) =-^ --^exp - / 

rit)) rit) \ do Hg 


dr 


zjr) 

r(T) 


exp dp^ dr + ^^^(O)^ . 


Note that both add and ^dd are bounded from below by Ki due to ()5.8I1 . We can then conclude that 7” is a 
self-mapping on 

M = Mo = L°“(0,r;L°“(G+)), 

and on 

M = Ml = L^iO, T- H\G+)) n L“(0, T; L“(G+)), 
since z G ^^(o, T; idi(G+)) n L°°(0, T; L°°(G+)). Moreover, 


{Gdw) — ^ ^^dd “t“ |d^(Zlid:,..., Dd—ldj^ ^) I ^ Pdd^ Iv—D^ih 


n' 


is small if || V'Ui||i°o(o,T;L“(ni)) and therefore || Vw||icxj(Q is small. This implies that T is a con¬ 

traction on 

Mil = {v G Mo ■■ 111) - Ddw||i«.(o,T;L“(G+)) < 7} 

and on 

Ml = {u G Ml : ||u - DdWIILcx,(o, t;L“(g+)) < 7} • 


13 
















for 7 sufficiently small. Thus the fixed point equation v = T{v) has a unique solution vg in Mg and it also 
has a uniqe solution vi in Mi , and both have to coincide vi = vg by uniqueness on Mg D Mi. On the other 
hand, obviously DdW lies in Mg, solves this fixed point equation and thus has to coincide with vg, hence 
also with vi- This proves that D^w € Mi C L^(0, T; 77^(G'+)). 

By transforming w back to u, we can conclude that u G The assertion then follows 

from the fact that the boundary is compact and can be covered by a finite set of balls {Br./ 2 {xi)}fLi. □ 

Higher boundary regularity was obtained under the assumption that u belongs to T; 

this was necessitated by the presence of the g-Laplace damping term in the equation. The assumption is 
equivalent to assuming Lipschitz continuity of u in time and space, i.e. u G T; (see Theorem 

4, Chapter 5, 0). 


5.1. Neumann problem for the coupled system 

It remains to show iJ^-regularity up to the boundary for the Neumann problem (113. 

Theorem 3. Let the assumptions of Corollary hold, let dfl and T = dLl.^. be regular and let u be 
the weak solution of (|4.5I) . Furthermore, assume that g G If u G T; IT^’°°(n)) 

and II V'ii||Lcx,(f 2 )(o_T;L°°(n)) sufficiently small, then Ui G H^{Lli)), i G {+,—}• 

Proof. We will show that u- G H^(0, T; since the regularity on 11+ follows as in the proof of 

Theorem [2 We begin again as before, by straightening the boundary around a;o G dfl- \i911+. There exists 
a ball B — Br{xg) for some r > 0 and a C^’^-diffeomorphism 'k : H ^ ^{B) C such that det|V'I'| = 1, 
[/' = n 11_) C M+ is an open set, and n 911_) C 91R+. We change the variables and write 

y = 4'(a:), xG B, 

X = y G U'. 

We denote G := '^{B^ n 11_) CC U'. We define 

w{y, f) := u($(y), 1), {y, t) G U' x [0, T], 

and transform the orginial equation from (B Oil) x [0, T] to U' x [0, T]: 


1 ^ 

(-{1 - 2kw{t))w{t)(j) + '^(aijDiw{t)Dj(j) + ^ijDiWi{t)Dj(j) 
lu' '' X . . W 


*j'=i 


|j|’Vw(t)|« ^ri,jDiw{t)Dj(j)) -^{w{t)f(l)\ dy = [ #| J,j, ^n|Rd dx, 

^ A ^ Jdu' 


(5.11) 


for a.e. t G [0,T], and all f G and A, k, u^, ^ij, and fjij are defined as in (15.21) . We then 

again use </> = —D~^{CfD\.w{f)), r G {1,..., d — 1} as a test function and proceed with the estimates like 
in the proof of Theorem [3 The only difference here is the need to estimate the boundary integral over dll' 
appearing in the weak form, therefore we focus our attention solely on estimating this term: 

— f f gDf'‘{('^Dlw)\J^^n\dxds 
Jo JdU' 

gDf\C‘^Dlu) dxds 


0 JBndn- 

< C'll5llL2(o,T;ffi/2(an_))IIC^^r^llL2(o,T;fl'i/2(an_)) 

< C'll5llL2(o,T;ffi/2(an_))||C-0[M||L2(o_r;//i(n_)) 

< e||Cdl(,'u||+ 2 (o,T;ffi(n)) + ^^^ll5llL2(0.T;ffi/2(ao_))- 

Showing that Dmw G H^{0,T; L‘^{G)) follows as in the proof of Theorem [3 since we use (f G G^(G) in 

(imi) . □ 
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